Berry Phase and Fidelity in the Dicke model with A 2 term 
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The instability, so-called the quantum-phase-like transition, in the Dicke model with a rotating- 
wave approximation for finite N atoms is investigated in terms of the Berry phase and the fidelity. 
It can be marked by the discontinuous behavior of these quantities as a function of the atom-field 
coupling parameter. Involving an additional field A 2 term, it is observed that the instability is not 
eliminated beyond the characteristic atom-field coupling parameter even for strong interaction of 
the bosonic fields, contrarily to the previous studies. 

PACS numbers: 42.50.Nn, 64.70.Tg, 03.65.Ud 



I. INTRODUCTION 

The Dicke model (DM) 1] describes an ensemble artifi- 
cial two-level atoms coupling with a cavity device. It has 
been attracted considerable attentions recently, mainly 
due to the fact that the Dicke model is closely related to 
many recent interesting fields in quantum optics and con- 
densed matter physics, such as the superradiant behavior 
by an ensemble of quantum dots Q and Bose-Einstein 
condensates [H , coupled arrays of optical cavities used to 
simulate and study the behavior of strongly correlated 
systems 0, and superconducting charge qubits[E [||. It 
is known from the previous studies 0,11, H that the full 
DM undergoes the second-order quantum phase transi- 
tion [nj. 

As claimed in Ref. [ll|, a sequence of instabilities, 
so-called quantum-phase-like transitions, is involved in 
the problem of an ensemble of two-level atoms system 
interacting with a bosonic field in the rotating-wave ap- 
proximation (RWA) [H, EE Ell • As addressed Ref. [l|, 
the absence of field A 2 term from the minimal coupling 
Hamiltonian in the approximation of the DM leads to 
the possibility of the instability. In the presence of A 2 
term, the classical thermodynamic properties have been 
studied previously [IE EE EH • Whether the instabilities 
disappear when the interaction of the bosonic field A 2 
term is taken into account is a long-standing issue and 
remains very controversial to datepE, EE Ell- 

It is known that quantum critical phenomena exhibits 
deep relations to the Berry phase _(BP) El EjU [HI E3 
and the fidelity [H, H, M, HE M, M, M, Hgj. The 
BP has been extensively studied by the geometric time 
evolution of a quantum system, providing means to de- 
tect the quantum effects and critical behavior, such as 
quantum jumps and collapse [IE HE [33[. A recent pro- 
posal is to use the fidelity in identifying the quantum 
phase transition [IE [H|. As a consequence of the dra- 
matic changes in the structure of the ground states, the 
fidelity should drop at critical points. In our opinion, the 
quantum-phase-like transitions might also be studied in 



terms of the BP and the GS fidelity. 

In this paper, we calculate the ground state BP and 
fidelity in the RWA DM with and without A 2 term to 
quantify phenomena of the quantum-phase-like transi- 
tions in finite system. Without A 2 , a exact solution to 
the DM is given explicitly. We solve the RWA DM with 
an additional A 2 term by a exact diagonalization in the 
Fock space of the bosonic operators. The paper is or- 
ganized as follows. In Sec. II, we review the RWA DM 
and the Hamiltonian with the A 2 term to obtain exact 
solutions respectively. In Sec. Ill, we study the instability 
of the RWA DM by measuring the BP and the ground 
state fidelity. The behaviors of these two quantities as 
a function of the interaction strength of the field for the 
RWA DM with A 2 term are also evaluated. Finally, we 
present the conclusion in Sec. IV. 



II. MODEL 
A. Exact solution to the RWA DM 

Let us consider DM of N two-level atoms with energy 
level wo j interacting with a single-mode bosonic field with 
the frequency u>. In the RWA DM, ignoring the counter- 
rotating term, the corresponding Hamiltonian is given 
by 



H = loo) a + luq,J z H -={a) J- + aj^ 

v N 



(1) 



where a) and a are the photonic creation and annihilation 
operators, Jk{k = z,±) denotes the collective spin-1/2 
atomic operators, A is the atom-field coupling strength, 
and h is set unity. 

Motivated by the exact technique of the Jaynes- 
Cammings model with RWA, we present a detailed nu- 
merical diagonalization procedures to solve a set of closed 
equations to the DM with RWA, which was also briefly 
discussed in Ref. .12}]. Since the Hamiltonian ( [T]) 



2 



commutes with the total excitation number operator 
L = a^a + J z + |, the subspace of the Hilbert space con- 
sists of a sum of subspaces labeled by different number of 
excitations L. The Hilbert space of the collective algebra 
is spanned by the kets to); to = —j, — j + 1, • • • , j}. 
By adapting Schwinger's re pre sentation of spin in terms 
of harmonic oscillators (34l. l35t| . \j,m) can be expressed 
as \j, N — n), which is a Dike state of N — n spin-up 
atoms and n spin-down atoms, n — 0, 1, ...N. In this 
work, j takes its maximal value N/2. \N/2,N — n) 
is also known as the eigenstates of J z and J 2 with 
J z \§,N-n) = (f -n)|f ,N-n). The action of the cor- 
responding raising and lowering operators on this state 
gives 



J+\—,N-n) = y/(N-n + l)n\ — ,N-n+l) 



J_\-,N-n) = y/(N-n)(n+l)\-,N-n-l) 
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FIG. 1: Excitation number L versus atom-field coupling pa- 
rameter A for different numbers of atoms N — 1,2, 4, 8. 



In the subspace of L 
function is supposed as 



N + k excitations the wave 



JY 



= '£c n \n + k) f (g)\N/2,N-n) a 

71 = 

(k = -N,-N+l,...) 



(2) 



where c„'s are coefficients, \n + k) f is a Fock state of the 
bosonic field with an alterative number k, ranging from 
—N to infinity, k is equal to — N in the weak coupling 
regions corresponding to excitations and then increases 
with the coupling parameter A. When the number of 
excitations L is larger than the number of atoms TV, i.e. 
k > 0, the ground state of H lies in the subspace spanned 
by N + 1 vectors. In this way, the Hamiltonian is ex- 
pressed by a tridiagonal (N + k+1) x (N + k+1) matrix. 
From Eq. ( [5]) we obtain the exact expression of the m-th 
row of the Schrodinger equation 



ECr, 



A 



Vto + k \J (N — m + l)mc m -i 



+ —== Vm + k + \yJ{N - m)(m + l)c m+1 
v N 



where 



+ [w(to + k) + uj a (N/2 - m)]c n 



0,1, ...N+k : fc<0 
0,1,...N : fc>0 



(3) 



(4) 



Note that the above equations are closed and the set 
of linear equations for c's takes a tridiagonal form. Solu- 
tions for a given k are readily obtained through Gaussian 
elimination and back substitution. Finally the chosen k 
corresponds to the lowest energy among eigenvalues of 
the solutions for a fixed coupling strength A. It is in- 
teresting to find that the excitation number L is added 



step by step and keeps a constant in a coupling parame- 
ter interval [Xi, Xj], where Xi(Xj) is a quantum-phase-like 
transition point, as shown in Fig. [T] The first transi- 
tion point is denoted as A°. The sensitive quantities like 
the ground state BP and the fidelity will be calculates to 
quantify the discontinuities, so called instability, in the 
finite DM with RWA. 



B. Numerical exact diagonalization to the RWA 
DM with A 2 term 

It is interesting to discuss the effect of the interacting 
bosonic field in the atom-field system. As the interac- 
tions vector potential A, caused by the longitudinal part 
of the bosonic field, are taken into account, the Hamil- 
tonian in the RWA DM can be evaluated with an addi- 
tional term A 2 [I^,[l3|.The additional term A 2 has been 
discussed classically about thermodynamic properties by 
Rzazewski et al [l7l . [l8j . To extensively quantify the 
contributions of the A 2 term, we employ the quantum 
information tools such as the BP and the ground state 
fidelity to detect the quantum-phase-like transitions. 

In terms of the bosonic operators, the A 2 term is given 
by e(a 1f + a) 2 , where e is the interacting strength of the 
bosonic field. The overall Hamiltonian of the ensemble 
of two-level atoms interacting with the bosonic field is 
expressed as 

H A = uja t a + uj J z + ^=(a t J- +a.J + )]+e(a t + a) 2 (5) 



In order to obtain the numerical exact solution, we per- 
form a standard Bogoliubov transformation by introduc- 
ing bosonic annihilation (creation) operator 6(6^), such 



that fet 



a, a 



and 



Ml' 



1. After substituting 



into Eq. ( [S]) the total Hamiltonian is diagonalized 
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H A = \/uj 2 + 4wetfb + lu J z + TjiV^ 2 + Aloe - 



ui) 



+ bJ + )-v{tfj + + bJ_)\ 



(6) 



where 
2 

M : 



2c 



V^ 5 



4lje 



-( — 



2s 



4we 



1) 



Note that a counter-rotating term is included in the mod- 
ified Hamiltonian ( [6]) , which may plays an essential role 
in the following discussion. Because the A 2 term breaks 
the gauge invariance of the Hamiltonian ( [T|) in the DM 
with RWA, it was argued in Ref. the instability 

would be then eliminated. 

We now consider the wave functions of the total Hamil- 
tonian with N atoms, which are of the form 



W)a = 



N Ntr 

^ ^ ^ ^ dnm 

n— m— 



m)f\N/2,n) a 



(7) 



where Ntr is the maximum photonic number in the ar- 
tificially truncated Fock space, and d nm are coefficients. 
\m) f is a Fock state with m photons. \N/2, n) a is a Dickc 
state in Schwinger's representation of spin with n atoms 
in excited state. The m-th row of the Schrodinger equa- 
tion reads 



N 1 

Ed nm = [uj £ m + A(n - —) + -{w e 



■4^=\/(m - 



1)(A- 



l)nd 



n — l,m+l 



-j= y/m(n + 1)(N - njd n +i,r, 
v N 



-j= \J{m + + 1)(A - n)d n+1 , m+ i 
v N 



Xv 



v / m(iV - 



l)nd r , 



(8) 



The eigenvalues and eigenfunctions can be obtained nu- 
merically by diagonalizing a (N + 1) x (Ntr + 1) matrix. 
The BP and the fidelity can be calculated through these 
eigenfunctions. 

To be complete, we also briefly review the contribution 
of the A 2 term in the DM model without RWA, which 
yields some unimportant corrections. The Hamiltonian 
of the full DM with the A 2 term reads 

Hdm = ojcfia + ujq,J z H — j=(a^ + a,)J x + e(a i + a) 2 (9) 



N 

With a rotation around an y axis by an angle \ and the 
same Bogoliubov transformation, the modified Hamilto- 
nian Hdm is rewritten as 

\J uj 2 + Aojetfb - uj J x + \^(\J J 2 + 4we - to) 



H 



DM 



2A 



(n-v)(tf+b)J, 



(10) 




FIG. 2: The average Berry phase jt/N of the RWA DM as 
a function of the coupling parameter A for different numbers 
of atoms N = 1,2,4,64. The inset shows a discontinuous 
picture of instability for 64 atoms. 



Therefore, A 2 term in the DM model without RWA would 
not change the nature of the phase transition, except that 
the position of the critical point is shifted. 



III. GROUND STATE PROPERTY 

A. Instability in the RWA DM 

Berry's adiabatic geometric phase describes a phase 
factor of the wave functions in a time-dependent quantum 
system. The interesting paths of evolution for generating 
a BP are those for which the ground state of the system 
can evolve around a region of criticality. We first mea- 
sure a nontrivial BP circulating a region of " criticality" 
corresponding to a abrupt change. The BP 71 generated 
after the system undergoing the time-dependent unitary 
transformation U(T) — exp[— i<f>(t) J z ], varying the angle 
cj)(t) adiabatically from to 2ir, can be evaluated as a 
function of the atom-field coupling parameter A 

7i - i r (MU\t)^-U{t)\^)d4> = 2^(Vo| JzWo) 
Jo d< P 

(11) 

where \tpo) is the ground-state wave function of Hamil- 
tonian m of the RWA DM. As shown in Fig. d the 
first phase-transition-like occurs at the "critical" value of 
the coupling parameter A" = 1 for arbitrary atom num- 
ber N, which recover the result in the thermodynamical 
limit [3a|. The average BP 71/iV is equal to n at A < 1 
and increases abruptly at discontinuous "critical" points 
when A > 1. Note that the plateau is formed clearly for 
N = 1,2,4 , the width of the plateau becomes narrower 
and narrower with the increasing A, which are quite dif- 
ferent from the phenomenon of the quantum phase tran- 
sition in the full DM @, H, H|. A clear picture of the 
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FIG. 3: The average Berry phase 72/iV.i.e. the average 
phonon number, in units of 2n of the RWA DM as a function 
of the coupling parameter A for different number of atoms N. 
The inset gives a discontinuous picture of instability for 64 
atoms. 



FIG. 4: Ground state fidelity F and energy gap A between the 
energies of the first excited and ground state in the RWA DM 
as a function of the coupling parameter of A/A° for different 
number of atoms N — 1,2, 4, 00. 



instability in the ground state of the RWA DM is given 
in terms of the BP with N — 64 atoms shown in the inset 
of Fig. El 

The effect of decoherence of the driving field on adia- 
batic evolutions of spin and quantized modes system has 
been investigated [lo, HH • In the fully quantized context 
we need a procedure capable of generating an analogous 
phase change in the state of the field. The BP 72 is ob- 
tained in terms of the bosonic operator by the phase shift 
unitary operator R(cf>) = exp[—i(j>(t)n], where h = a^a is 
the number of bosons in the field. Changing the angle 
4>{t) slowly from to 2ir the ground state 72 is given by 

72 = i [ {^ \R\t)^-R{t)\^)dct> = 27r(^ | a t a |^ ) 
Jo # 

We now have a general expression for the BP 72 re- 
lated to the photonic number, which is driven by fields. 
We plot behaviors of 72 /N in units of 2ir as a function 
of the atom-field coupling parameter A for different num- 
ber of atoms N in Fig. \3\ The BP j 2 /N is in the 
weak coupling region for A < 1 and then increases dis- 
continuously as A increases. As shown in Fig. [3J when 
N increases the interval of the "critical" values A become 
smaller, leading to the curve with more steps. The inset 
of Fig. [3J shows that there actually exist many phase- 
transition-like "critical" points beyond A = 1 for large 
N = 64. 

An increasing interest has been drawn in the role of 
the ground state fidelity in detecting the quantum phase 
transitions for various many-body system, with a narrow 
drop at the transition point. Below we propose to use this 
quantum tool to identify the quantum-phase-like transi- 
tions, where the GS fidelity drops to in the RWA DM. 
It is defined as the overlap between two ground states 
|^o(A)) and ^(A + fJA)) [2g, |27|, where <5A is a tiny per- 



turbation parameter, i.e. 

F(\,5\) = \(M*)\MX + 5\))\ (13) 

Note that F is a function of both A and <5A. Based on the 
normalized and orthogonalized wave function in Eq.([2]) 
for the RWA DM, the ground state fidelity can be derived 
analytically | J2n, m =o c„(A)c m (A + 5\)5 n>m 5 n+kim+k < |, 
and then can be simplified as 

F(A,<5A) = (° : n = m,k^k (M) 
L_ 1 : n = m,k = k 

At each transition point in RWA DM, the alternative 
number k is changed abruptly and F is then equal to 0. 
Beyond "critical" points, F should be constant 1. The 
energy gap A between the first excited and ground state 
energies tends to at each "critical" point, as shown 
in Fig. [U We attribute this level crossing to the fact 
that the GS wave functions at the different sides of each 
transition point arc orthogonal. The ground state fidelity 
drops to at each "critical" point, exhibiting similar crit- 
ical singularity of the first-order quantum phase transi- 
tion [2g|. It is also obvious that the instabilities increase 
with a increasing number of atoms. 



B. Behaviors of the Berry phase and fidelity in the 
RWA DM with A 2 term 

We next turn to study the RWA DM with A 2 term by 
using the above quantum information tools. By manse of 
the general BP formula Eqs.f [TTj) and ( IT!?]), the ground 
state BP of the RWA DM with A 2 term can be evalu- 
ated as a function of the atom-field coupling A and the 
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FIG. 5: The average Berry phase (a) jia/N, (b) ')2a/N in 
units of 2-7T, as well as (c) ground state fidelity Fa of the RWA 
DM with A 2 term versus A for different interacting strengthes 
of the field e = 0, 0.5, 1 for N = 4. 

interaction strength of the field e, i. e., 

Iia = 2n((p A \J z \(p A ), 

l2A = 2n(ip A \bh\ip A ) (15) 

where \(p A ) is the ground state of the Hamiltonian (J5]). 

We plot the average BP ji A /N in Fig. [5](a) and ^a/N 
in units of 2it in Fig. 0(b) as a function of the coupling A 
for different interaction strengthes e = 0, 0.5, 1 for N = 4 
atoms. We can observe that, for both ji A /N and ^za/N , 
the abrupt jump occurs at the same coupling parameter 
A for the same interaction strengthes e. We denote the 
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FIG. 6: Ground state fidelity Fa (a), Berry phase Jia/N (b) 
and its first derivative d'yiA/dX (c) of the RWA DM with 
A 2 term as a function of the coupling of A/A c for different 
interacting strengthes of the field e = 0, 10 -3 , 10 -2 , 10 -1 , 1. 



coupling parameter A where the first abrupt jump oc- 
curs as a characteristic parameter A c . It is interesting 
to observe that the characteristic A c increases with the 
interaction strength e. The stability claimed in Ref. [l5[ 
is only found in the smooth curves for A < A c . Thus, 
there still exist quantum-phase-like transitions when the 
interaction of the bosonic field e is strong. 

According to the ground state wave function \fA), the 
fidelity of the RWA DM with A 2 term can be also calcu- 
lated ,Fa(A,<5A) = \((p A (\)\(f A (X + S\))\. The numerical 
results for the different e with N = 4 atoms are exhibited 
in Fig. [5Jc) . The singularities at the " critical" points for 
e = 0, 0.5, 1 are demonstrated by a sudden drop of F A . 
It is clearly shown that the characteristic A c moves to- 
wards the right regime with the increasing e, providing 
the evidence of the quantum-phase-like transitions even 
for a strong interaction of the field. 

To show the instabilities more obviously, for a 
more wide range of interacting strengthes e = 
0, 10~ 3 , 10 -2 , 10 _1 , 1 , we plot the ground state fidelity 
F A , Berry phase ~f\ A /N and its first derivative dji A /d\ 
of the RWA DM with A 2 term as a function of the scaled 
coupling parameter A/A c together in Fig. [BJ When the 
interaction strength of the field increases the fidelity F 
still drops to at the characteristic A c , where the first 
derivative N~ 1 dji A /dX also changes abruptly. This is 
another piece of evidence that the instability of the RWA 
DM does not vanish in the RWA DM including the A 2 
term. 

It is illustrated that the contribution of the A 2 term 
does not eliminates the instability of the RWA DM, con- 
trarily to the previous studies by Rzazewski et al.[ll|. 
For strong interaction strength of the bosonic field e, a 
characteristic parameter A c becomes larger than A° = 1 
in the absence of A 2 term. A sequence of the ground state 
stability reported previously only appears for A < A c . 
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IV. CONCLUSIONS 

In summary, we have investigated the instability of the 
RWA DM by quantum information tools such as the BP 
and the ground state fidelity. An obvious discontinuous 
behaviors of these quantities with finite N atoms are ob- 
served. It is demonstrated that the quantum-phasc-likc 
transitions occur beyond the characteristic A c for strong 
interaction of fields. We propose that the instability 
would not be eliminated by involving the A 2 term of the 
DM with RWA. Previous observed instability may be lim- 
ited to the coupling regime A < A c . It should be pointed 
out that the quantum information tools are very sensitive 
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